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Spin-up of low mass classical bulges in barred galaxies 
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ABSTRACT 

Secular evolution is one of the key routes through which galaxies evolve along the 
Hubble sequence. Not only the disk undergoes morphological and kinematic changes, 
but also a preexisting classical bulge may be dynamically changed by the secular 
processes driven primarily by the bar. We study the influence of a growing bar on the 
dynamical evolution of a low mass classical bulge such as might be present in galaxies 
like the Milky Way. Using self-consistent high resolution A^-body simulations, we study 
how an initially isotropic non-rotating small classical bulge absorbs angular momentum 
emitted by the bar. The basic mechanism of this angular momentum exchange is 
through resonances and a considerable fraction of the angular momentum is channeled 
through Lagrange point (-1:1) and ILR (2:1) orbits. In the phase of rapid dynamical 
growth, also retrograde non-resonant orbits absorb significant angular momentum. 
As a result of this angular momentum gain, the initially non-rotating classical bulge 
transforms into a fast rotating, radially anisotropic and triaxial object, embedded in 
the similarly fast rotating boxy bulge formed from the disk. Towards the end of the 
evolution, the classical bulge develops cylindrical rotation. By that time, its inner 
regions host a "classical bulge-bar" whose distinct kinematics could serve as direct 
observational evidence for the secular evolution in the galaxy. Some implications of 
these results are discussed briefly. 
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- galaxies: spiral - galaxies: evolution 



1 INTRODUCTION 

In the hier a rchica l structure formation scenario 
(|White fc Reesl 1 19781 : iFall fc Efstathioul Il980h mergers 
have played a strong role in forming and shaping galaxies. 
One of the co mmon product of major mergers are the 
classical bulges (iKauffmann et all Il993l : iBaugh et all 1 19961 : 
iHopk ins et al.ll2009l ) which are the central building blocks 
in spiral galaxies. There have been a couple of other 
mechanisms suggested for the formation and growth of 
classical bu lges e.g., monolit hic collapse of primordial 
gas clouds (|Eggen et al.l 1 19621 ). the coalesc ence of giant 
clum p s in gas-rich primord ial galaxies l|lmmeli et al.l 



sense as the disk rotates. It is also known that classical 
bulges rotate faster than elliptical galaxies and that often 
their rotation velocities are comparable t o that of an 
isotropic oblate rotator model (|Binnevlll978f ). So the origin 
of the systematic rotational motion observed in the classical 
bulges remains unclear. 

The photometric and kinematic properties of classical 
bulges as well as their origin are quite distinct from those 
of the other class of bulges, the boxy/peanut and disk- like 
bulges. It is well known that the surface brightness profiles 



with sersic 
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.20081). multiple m inor mergers 
i Bournaud et al.1 12007^ iHopkins et al.l |2010|), accretion of 
small companions or satellites (|Aguerri et all 120011 ). The 
classical bulges formed via these processes seem to have 
little rotation as compared to the random motion. On 
the other hand, various observational measurements have 
confirmed that classical bulges in spira l galaxies possess 
rotat ion (|Kormendv fc Illingworthl 1 19821 : ICappellari et al] 
|j'H)7l i about their minor axis and in most cases in the same 
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in classical bulges follow a sersic law p(r) 
index n ~ 4. Whereas the sersic indices in boxy and disk-like 
bulges are, in general, low with n ^ 2; so that their surface 
brightne ss profiles follow roughly an e x ponentia l distr ibu- 
tion, see lKormendv fc Kennicutrj |2004h : bombed (|2009l ) for 
extensive reviews. The kine matics of bulges are well illus- 
trated in the v/ a — e plot (|Kormendv fc Illingworthl Il982l ; 
lKormendvlll982n which clearly demonstrates the distinction 
in the kinematic properties of ellipticals, classical bulges and 
boxy bulges and brings out the fact that in terms of their 
rotational support, classical bulges fall in between ellipticals 
and boxy /disk-like bulges. Since the boxy as well as disk- like 
bulges are thought to have formed from disk material, the 
source of their angular momentum is known. 
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The classical bulges formed early through major merg- 
ers and violent relaxation can subsequently accrete mate- 
rial quiescently as a result of which a disk grows inside- 
out Mo et al. Hl998t iKatz et alJl2003l : ISpringel fc Hernquistl 
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Keres et al. 20091 ). The gas accretion may facilitate the 



disk to grow sufficiently for the disk self-gravity to dominate 
the internal dynamics. Eventually, a bar and/or spiral arms 
form in the disk and initiate secular processes in the galaxy. 
Indeed, bars are quite common in disk galaxies, about 2/3 
of the disk galaxies host a strong stellar bar in their cen- 
tral region (lLaurikainen et al.ll2004l ; lMarinova fc Jogee|[2007l ; 
iMenendez-Delmestre et al.ll2007f K Therefore one might ex- 
pect the disk to go bar unstable also in galaxies with pre- 
existing classical bulges. I n fact, classical, boxy and disk- 
like bulges could co-exist (lAthanassoulal 120051 ; lErwmll20Qg| ; 
iGadottil 120091 ; iNowak et all |2010T | in a single galaxy, al- 
though the observational identification of the several com- 
ponents could be difficult. In the Milky Way, an upper limit 
on the mass of a classical bulge (~ 8% of the disk mass) has 
been set by modelling the kinemati cs from the Bulg e Ra- 
dial Velocity Assay (BRAVA) data (|Shen et alj|2010l ). But 
there is ev idence for a metallicity gradient a bove the Galac- 
tic plane (|Zoccali et al.l [200 i lZoccali|[2010l ) which is taken 
as an indication for the existence of a classical bulge in our 
Galaxy. It is therefore important to understand the dynami- 
cal interaction between preexisting classical bulges and bars 
in barred galaxies. 

In this paper, we investigate in considerable detail the 
interaction of a bar and a low mass classical bulge via a high 
resolution TV-body simulation of a galaxy consisting of a live 
disk, bulge and dark matter halo, and follow the evolution of 
the dynamical structure and kinematics of the small classical 
bulge. We find that its dynamical evolution is strongly con- 
nected to the growth of the bar which forms spontaneously 
in the disk. During the secular evolution, the structure and 
kinematics of the bulge are altered significantly, developing 
an interesting and complex rotation structure; in particular, 
cylindrical rotation (which is considered as a typical proxy 
of boxy bulge) appears in the inner region of the classical 
bulge. 

The paper is organized as follows. Section[2]summarizes 
the basics of bar-bulge interaction. Section[3]outlines the ini- 
tial galaxy model and set up for the TV-body simulation. The 
bar evolution and boxy bulge formation are described in Sec- 
tionfil Section[S]describes, in detail, the angular momentum 
exchange between the bar and the classical bulge. The evo- 
lution of the classical bulge, its structure and kinematics are 
presented in Section [5] Discussion and conclusions are con- 
tained in Sections[7]and[5]respectively. In the text, by bulge, 
we mean a classical bulge unless mentioned otherwise. 



2 BAR-BULGE INTERACTION 

As we have seen (in Section [1} , the possible co-existence 
of a bar and a classical bulge might be rather common in 
present day disk galaxies, and thus they are bound to inter- 
act gravitationally. In fact, a preexisting classical bulge in 
the disk has a strong influence on the formation and growth 
of the bar itself. For the swing amplification to work, one 
needs to keep alive the feedback loop through which a set of 
trailing waves traveling through the center are transformed 



into leading waves. This is possible as long as there is no 
inner Lindblad resonances (ILRs). A highly centrally con- 
centrated bulge can shield the center by putting an ILR 
barrier and thus cutting the feed back loop which in tur n 
could hinder the growth of the bar (|Sellwood fc Evansl200ll '). 
However, various non-linear processes are probably active in 
real galaxies which would destroy the I LR barrier and even- 
tually lead to the fo rmation of a bar jWidrow et al.l [20081 ; 
iDubinski et al.ll2009T ). 

Once a bar is formed, it takes over the dynamics in the 
central region of the disk and starts interacting with the stel- 
lar bulge and dark matter halo thro ugh exchange of angular 
mome ntum. Based on the work of iLvnden-Bell fc Kalnaia 
(| 19721 ) , he reafter LBK72, it has b e en em p hasized by s evera l 
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resonant interaction plays a significant role in the angular 
momentum transfer between the b ar and the dark matte r 
halo. It has been suggested by iHernquist fc Weinber 

[Weinberg fc Katzl fcOOZ 



It has been s u ggeste d 
(|l992T ); lAthanassoulal (20031 ); 



I) 

that the same underlying mechanism could as well apply 
between the bar and th e sph eroid and in particular, 
lAthanassoula fc Misiriotisl |2002T ) have studied how the 
shape of a bulge would change in response to a growing bar. 

Although the dynamical interaction between a growing 
bar and a bulge and their subsequent evolution can be best 
studied via TV-body simulations, an analytic understanding 
is required to complement this. Following LBK72, it can be 
shown that during the bar-bulge interaction the time rate 
of change of angular momentum of a classical bulge, whose 
distribution function (Fb) is descibed by a King model, is 
always positive and can be written as 



L Z)h ~ x \i>i mn \ x F b /al > 0, 



(1) 



where ijjimn and Qb are the Fourier amplitude and pat- 
tern speed of a non-responsive bar potential and ctj, is the 
velocity dispersion of the bulge stars. So at a given reso- 
nance, the angular momentum gained by the bulge depends 
on the strength of the bar and is inversely proportional to 
the square of bulge velocity dispersion, implying a hotter 
bulge will absorb less angular momentum provided other 
conditions remain unchanged. However, in real galaxies, the 
angular momentum transfer is more difficult to determine, 
because the time rate of change of bar's angular momentum 
involves the change in its pattern speed, moment of inertia, 
and in the ang ular momentum associat ed with any inter- 
nal circulation (| Villa- Vargas et all 120091 ) within the bar. In 
Section [SJ we show the angular momentum transfer between 
the bar and the bulge in our simulation using orbital spectral 
analysis. 



3 GALAXY MODEL AND TV-BODY 
SIMULATION 

An equilibrium model for a disk galaxy is con- 
st ructed using the se l f-cons istent bulge-disk-halo model 
of iKuiiken fc Dubinskil (1 19951 ). Their prescription provides 
nearly exact solutions of the collisionless Boltzmann and 



Spin-up of low mass classical bulges 



Poisson equations which are suitable for studying disk sta- 
bility related problems. All the components in our model 
are live (i.e., the gravitational potential of each component 
can respond to an external or internal perturbation) and, 
hence, provide a realistic representation for the structure 
and evolution of the galaxy. Below we briefly describe each 
compone nt of the model. For more d etails, the reader is re- 
ferred to lKuiiken fc Dubinskil (| 19951 ). 

The disk distribution function is constructed using the 
approximate third integral given by E z = + ty(R,z) — 
ty(R,0), the energy of the vertical oscillations. This third 
integral is approximately conserved for orbits near the disk 
mid-plane. The radial density of the disk is approximately 
exponential with a truncation, and the square of the radial 
velocity dispersion follows the same exponential radial de- 
cline with a scale length same as the disc scale length. The 
vertical structure of the disk is approximately isothermal, 
with the scale height set by the vertical velocity dispersion 
and vertical potential gradient. The volume density of the 
axisymmetric disk is given by 



Pd(R,z) 



M d 



inh z R d 



-R/Rd 



erfc 



R Ron 



^fi{Rout — Rt 



M*),(2) 



where f d (z) = exp(-0.8676*z(.R, z)/V z (R, h x )) with 
$ z(R, z) — ^(R, z) — ^(R, 0) governs the vertical structure 
of the disk, erfc is the complementary error function. In the 
above equation, Md is the disk mass, Rd is the scale length 
and h z is the scale height. 

A spher ical live cla ssical bulge is constructed from the 
King model (|Kinelll96a ) a nd the corresponding distr ibution 
function (DF) is given by jBinnev fc Tremaine 1 19871 ^ 



ME) 



p 6 (27T<7 b 2 )- 3 /V*«>-*^ CT fc 
x { e -(B-*c)/<r= _ if£<* c , 

otherwise. 



(3) 



Here, the bulge is specified by three parameters, namely 
the cut-off potential which determines the bulge tidal 
radius), central bulge density (pb) and central bulge velocity 
dispersion (<7(,). The gravitational potential at the centre of 
the bulge is measured by ^bo- 

An axisymmetric live dark matter halo is constructed 
using the di stribution function of a lowered Evans model 
jEvansil 1993T ) and is given as 



fdm(E,L z ) = 




if E < 0, 
otherwise. 



(4) 



The halo is parameterized by a potential depth (^o), veloc- 
ity (an) and density scales (pi), a core radius R c and the 
flattening parameter q. The factors A, B, and C are fu nc- 
tions of these parameters fsee lKuiiken fc Dub inski 1993 an d 
references therein). The halo has a tidal radius specified by 
E = 0. 

The total mass and the outer radii of both the bulge 
and halo are calculated in an iterative procedure. The po- 
tential is computed self-consistently by solving the Poisson 
equation for the combined three component system in an 
iterative fashion. First, the densities for the bulge and halo 
are obtained from their respective distribution function and 
then the disk density is added to it and the corresponding 
potential for the combined mass distribution is used as a 
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Figure 1. Initial circular velocity curve for the model galaxy. 
Solid line represents the total circular velocity . Dotted line is for 
the bulge, dashed line for the dark matter halo and dashcd-dot- 
dash for the disc. 



starting point for carrying out the next iteration. We use 
a maximum of I — 10 in the potential harmonic expansion 
and the iteration is continued until the outer radii for the 
bulge and halo are unchanged between successive iterations. 
The outer radii of the bulge and the halo correspond to the 
respective tidal radii. The masses of the bulge and halo cor- 
respond to the total mass enclosed within their respective 
outer radii computed by integrating the density profiles. 

In this paper, we present the analysis of a particular 
galaxy model hosting a low mass classical bulge. For histori- 
cal reasons, we call this model RCG004. The circular velocity 
curve for the model is presented in Fig[T] The length, mass 
and velocity units for this model are given by L = 4.0 kpc, 
M = 2.33 x 1O 1O M and V = 157 kms" 1 . The disk outer 
radius (R ou t) is fixed at about 6.5Rd and a truncation width 
~ 0.3Rd is adopted within which the disk density smoothly 
decreases to zero at the outer radius. The disk scale length 
(Rd) is fixed at 4.0 kpc and the scale height is 42 pc, the disk 
mass Aid = 4.5 x 10 10 M Q . The central value of the radial 
velocity dispersion is 78.5 kms -1 . The Toomre Q profile is 
nearly flat in the radial range 0.5 to 5 scale lengths while it 
increases on either side of the disc. The Q value at the disc 
half mass radius is 1.4. The bulge mass Mb = 3 x 10 9 Mq. 
In table [TJ we quote the outer radius for the classical bulge 
(denoted by Rb) in our galaxy model. The halo has a flat- 
tening of q — 0.8 and a core radius R c — 0.25 kpc and a 
mass of M h = 1.82 x IO^Mq within about 60 kpc. 

We evolve the galaxy model in isolation to examine 
the evolution of the bulge shape, morphology and kine- 
matics. The simulati on is performed using the Gadget code 
(|Springel et al.ll200ll ) which uses a variant of the leapfrog 
method for the time integration. The forces between the 
particles are calculated using the Barnes fc H ut (BH) tree 
with some modification (|Springel et al.l [200lh with a toler- 
ance parameter 9 to i = 0.7. The integration time step is ~ 0.4 
Myr and the model is evolved for 2.2 Gyr. For reference, the 
orbital time at the disk half mass radius is ~ 296 Myr. A to- 
tal of 1.0 x 10 7 particles is used to simulate the model galaxy. 
The softening lengths for the disk, bulge and halo particles 
are 12, 40 and 36 pc respectively. The masses for the disk, 
bulge and halo particles are 1.2 x lO 4 Af , 0.3 x 10 4 Af Q and 
3.6 x 1O 4 A40 respectively. To examine the effect of unequal 
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Figure 2. Surface density maps for the disk particles alone. Top 
left: density map at T=0 Gyr, showing the axisymmetric disk. 
Top right: Same at T=0.56 Gyr. Bottom left: at T=l.l Gyr and 
Bottom right: at T=2.1 Gyr. 



Table 1. Initial parameters for the model galaxy. 



Galaxy 


Q B/D 


B/T a bQ 


Rb 


model 




(kms- 1 ) 


(kpc) 


RCG004 


1.40 0.0666 


0.01306 65.0 


6.08 



B/D is the bulge-to-disk mass ratio, B/T is the bulge-to-total 
(including dark halo mass) mass ratio, tr^o is the bulge central 
velocity dispersion , and finally ijj, is the outer radius of the bulge. 



softenings, we have re-run the simulation with new soften- 
ing parameters as prescribed bv lMcMillan fe Dehnenl (|2007t ). 
We note that with the new softenings, the bar growth is de- 
layed by ~ 90 Myr while the main results remain unchanged. 
The total energy is conserved within 0.2% till the end of 
the simulation. The total angular momentum is conserved 
within 3% at 2.2 Gyr for both the runs having different soft- 
ening parameters. 



4 BAR AND BOXY BULGE 

Although it is not clearly understoo d how bars are formed 
in real galaxies, swing amplification (|Toomrelll98ll ) plays a 
significant role in making an initially axisy mmetric, equilib - 
rium model of a disk galaxy bar unstable (|Sellwoodlll98ll ). 
Once formed, N-hoAj bars are found to be long-lived, domi- 
nate the disk dynamics, and are resp onsible for driving secu- 
lar ev olution processes in the galaxy (|Sellwood fc Wilkinson! 
Il995t l. Fig. [2] depicts the formation and evolution of the bar 
from the initially axisymmetric disk. Strong two-armed spi- 
rals are also formed along with the bar and last till 1.1 Gyr 
in our simulation. The ring-like structure at T = 0.56 Gyr 
is intersecting the spiral arms indicating that it is probably 




t (Gyr) 

Figure 3. Time evolution of the bar amplitude and the pattern 
speed (f2fl). Red solid line is the result of a linear regression 
analysis done on the measured pattern speed values from the N- 
body snapshots. 



not real. Such a ring-like feature arises because of the galaxy 
model not being in perfect equilibrium. 

In the upper panel of Fig. [3] we show the time evolution 
of the bar amplitude measured as the maximum of m = 2 
Fourier coefficient (A%) of the density perturbation normal- 
ized to the unperturbed axisymmetric component (^4o)- The 
bar reaches its first peak in amplitude at 0.28 Gyr and the 
second peak at 0.44 Gyr. The m = 1 vertical Fourier mode 
in the r — z plane corotating with the bar pattern 
speed shows that the disk is undergoing buckling instability 
from ~ 0.39 — 0.6 Gyr and strong buckling occurs around 
0.6 Gyr. 

Based on the nature of the growth curve, bars are clas- 
sified into two broad categories, type-I and type-II. Type-I 
bars are strong, grow within a few orbital time scales and 
nearly saturate in amplitude, whereas type-II bars are weak, 
the growth time scale is very long (typically, a se cular evolu- 
tion time scale) and show no sign of saturation |Saha et al.l 
2010). The bar in our model is a type-I bar (e.g., Fig. [3J. 

The bottom panel of Fig. [3] shows the evolution of the 
bar pattern speed. The pattern speed of the bar decreases 
over time in o ur simulation, primarily because of the dynam - 
ical friction (|Tremaine fc Weinberg! 1 19841 : IWeinberel 1 19851 ) 
against the dark matter halo. A detailed account on the 
bar's pattern speed decrease and its dependence on vari- 
ous dark matter halo properties e.g. halo angular momen- 
tum, orbital anisotropy, central concentration can be found 
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in iDebattista fc Sellwoodl (|2000f ) . Using a linear regression 
analysis on the simulation data, we find the half- life, T1/2, 
(the time period over which the bar pattern speed would 
decay to half its initial value) of the rotating bar to be 
~ 3.09 Gyr. This indicates that the rate of angular momen- 
tum transfer from the bar is rather slow in our simulation; 
for an in-d epth analysis on the bar slow down , reader s are 
referred to I Weinberg! (| 19851 ); I Weinberg fc Kat3 (|2007al ). 

As the bar grows stronger, its self-gravity increases 
and it goes through the w e ll-known buckling instabil 
itv dCombes fc Sandersl Il98ll; IPfenniger fc Normanl Il99 ' 



iRaha et al.1 Il99ll ; iMartinez-Valpuesta fc Shlosmanl 120041 ') 
following which the bar transforms into a boxy/peanut 
bulge. In Fig. [4] we present the surface density (left pan- 
els) and velocity field (right panels) for the boxy bulge seen 
edge-on; i.e., only disk particles are shown. The cylindrical 
rotation is evident. The final boxy bulge contains approxi- 
mately 33% of the disk mass including the inner barred disk 
component. Note that the density drops off sharply along 
the vertical direction in the boxy bulge region. In Section [6j 
we will compare the structure and kinematics of the boxy 
bulge in Fig. [4] formed in our simulation with the classical 
bulge undergoing the bar driven secular evolution. 



o 



o 



< 
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Figure 5. Evolution of the specific angular momentum of the 
bulge (green), disk (red) and halo (blue) components in our 
model. Along the y-axis plotted are the specific angular momen- 
tum minus its value at T = normalized by the disk angular 
momentum (L [ j z (0)) at T = 0. 



5 ANGULAR MOMENTUM TRANSFER TO 
THE CLASSICAL BULGE 

We compute the specific angular momentum for each species 
e.g., disk, bulge and halo particles in our simulation and 
re-confirm the already established fact that the inner re- 
gions of the disk loose angular momentum through the bar. 
While a significant fraction of the total angular momentum 
emitted by the bar is absorbed by the surrounding dark 
matter halo, the angular momentum gained by the bulge 
is non-negligible. In Fig.JS] we show the angular momentum 
transfer amongst the disk, bulge and halo components in our 
model. The total angular momentum is conserved within 3% 
at the end of 2.2 Gyr in our simulation. Initially both the 
bulge and halo have zero net angular momentum i.e., they 
start as non-rotating objects. Note that the rate of gain of 



angular momentum by the classical bulge particles nearly 
saturates towards the end of the simulation and closely fol- 
lows the growth of the bar (see Fig. [3}- Using orbital spectral 
analysis, we show below that the gain of angular momentum 
by the bulge occurs primari ly through resonances (see also 
iHernquist fc Weinberg]|l992f ). 

The simulation presented here shows an increase in the 
bulge rotation velocity (Section 16. 2[) . and a corresponding 
increase in bulge angular momentum. The transfer of angu- 
lar momentum from the bar to the bulge depends strongly 
on the pattern speed which sets the resonance locations. It 
is important to note that if the resonances are sparsely pop- 
ulated because of lack of particles in the si mulation, the an- 
gular momentum transfer will be inefficient (|Weinberg||l985l ; 
I Weinberg fc Katzl l2007bh . In our case, we have a total of 
10 7 particles with 10 6 particles in the classical bulge. So we 
can test whether angular momentum transfer through reso- 
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Figure 6. The top panels show the distribution of particles with frequency (CI — at five different times throughout the evolution of 

the bulge. The lower panels show the gain of angular momentum of the selected particles with respect to the previous time, as indicated 
on the top of the panel. The vertical dotted lines in the figure indicate the most important resonances, -1:1, 2:1, 3:1 and 4:1. Note that 
from time 1.1 Gyr to 2.1 Gyr the particles trapped at around (f2 — nt,)/K=0.5 gain less angular momentum than earlier. 



nant interaction, is the mechanism for the angular momen- 
tum gain of the classical bulge. Here, we quantify its effect 
by using an orbital spectral an alysis method described in 



Mart inez- Valpuesta et al.l (|2006r i based on that presented in 
Binnev fc Spergell (|l982j ). In previous works, this method 



was applied to halo resonant orbits ( Athanassoulal 120031 ; 
iMartinez- Valpuesta et al.ll2006l ; Ibubinski et alj|2009h to un- 
derstand the bar-halo interaction. We apply it to the classi- 
cal bulge particles to find out how many of them are trapped 
in resonances and what is the corresponding gain of angu- 
lar momentum. The potential is extracted from the Af-body 
simulati on at different snap s hots using the grid code pro- 
vided bv lSellwood fc Valluril (|l997h and then frozen to com- 
pute the orbits. We randomly select 100000 (10%) particles 
out of 1 million in the bulge and compute their correspond- 
ing orbits. We calculate the azimuthal and radial epicyclic 
frequencies Q and k respectively for each of the orbits by 
Fourier analysis. 

We present the results of this orbital spectral analysis in 
Fig.[5]at 5 different epochs. On the top panels, we present the 
classification of classical bulge particles by their frequency 
ratio r\ = (f2 — f2s)//t- The bar has an irregular evolution 
and this can be seen in the time sequence of the top panel in 
Fig. [6] Initially, the bulge particles are distributed half co- 
rotating with the disk and half counter-rotating. Therefore, 
when we study the orbital distribution in the very early stage 
of the bar growth at T = 0.17 Gyr, there still is an almost 
symmetric distribution. When the bar is already formed, 
right after reaching the maximum, at T — 0.39 Gyr, many 
particles have been trapped around the 2:1 resonance. Tak- 
ing a careful look at these orbits, we have checked that they 
are of xl-type. There is also a considerable group of particles 
with 77 £ (—0.4,0.); a look to the orbits allows to identify 
them as mainly stochastic trajectories. In the lower panels, 



we show the angular momentum gain by the particles at 
each 77 during the growth and evolution of the bar. There 
is a considerable gain of angular momentum by three main 
groups in our diagram (Fig. [6] second bottom panel). The 
main gain of angular momentum comes from those parti- 
cles at resonance with r] = — 1, corresponding to particles 
orbiting around the Lagrangian points. Since these particles 
are at negative frequency it means that they are counter- 
rotating with the bar. Another gaining group corresponds 
to the particles with r\ € (—0.4,0.), the stochastic group. 
By gaining angular momentum their (counter) rotation de- 
creases. Amongst the low order resonances, the important 
gaining group is around the ILR (77 = 0.5). We can conclude 
that at this stage of evolution, which is very rapid, the main 
transfer of angular momentum occurs through resonant and 
stochastic orbits. 

During the period of evolution between T = 0.56 — 
0.39 Gyr, the bar goes through the buckling event, there- 
fore there is still some trapping of particles around the 2:1 
resonance. Although there not much gain or loss of angu- 
lar momentum, some angular momenta are gained through 
r) = — 1 resonance during this period. Notice that in the 
upper row of panels, the number of bulge particles trapped 
at the bar's ILR (2:1) is increasing. In the forth panel, at 
T = 1.1 Gyr, we have 36% of particles trapped around the 
2:1 resonance. At this time, the main gain of angular mo- 
mentum comes from resonances at the ILR (rj = 0.5), and 
at 77 = —1,-2. The particles at the OLR (77 = —0.5) and 
those with 77 6 (—0.4,0.) are losing angular momentum. At 
T = 2.1 Gyr the classical bulge is still gaining angular mo- 
mentum through resonances at ILR and 77 = — 1 correspond- 
ing to the Lagrange points. Note that some of the angular 
momentum gain is also coming from the OLR. Although the 
number of particles trapped at the ILR is gradually increas- 
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ing over time, their angular momentum gain does not follow 
accordingly. By comparing the last two panels (upper and 
lower), it is evident that the particles trapped at the ILR 
are now hardly gaining angular momentum. It is plausible 
that the inner bar-like structure in the classical bulge (see 
section 153)) is giving away angular momentum to the outer 
parts of bulge and perhaps to disk and halo. 

The gain of angular momentum by the bulge can thus 
be explained by resonances during the slow secular evolu- 
tion of the bar, and by resonances together with stochastic 
orbits in the dynamical stage. During the dynamical phase, 
T — 0.56 — 0.17 Gyr, the net gain of angular momentum 
(computed by adding up the averaged angular momentum 
of each orbit) is 3 times larger than that gained in the rela- 
tively quiet secular phase (T = 2.1 — 0.56 Gyr). While ap- 
proximately 70% of the net angular momentum gain comes 
from the resonances, stochastic orbits contribute to ~ 30% 
of the net angular momentum gained during the dynamical 
phase. 

Previous studies of angular momentum tra n sfer t o 
the live dark matter h alo (e.g., Athanassoula (2003); 



iMartinez-Valpuesta et ail (|2006h : iDubinski etafl (120091 ) ) 
have found important contribution from corotation, outer 
Lindblad resonance (OLR) and higher order resonances. By 
contrast, in the case of a small classical bulge as studied here, 
the OLR and corotation have not played any significant role 



in the gain of angular momentum as shown above. This is 
most probably due to the fact that the size of the bulge in our 
simulation is much smaller than the typical size of the dark 
matter halo; the bulge half-mass radius (Ri/ 2 = 0.21_R<j and 
0.225i?d at T — 1.1 and 2.1 Gyr respectively) is shorter than 
the bar size (Rbar = 0.987i?,j and 2.01-Rd at those times) in 
our simulation. The bar size is measured from t he phase an- 
gle of the bar (I Athanassoula &i Misiriotisll2002r i. The phase 
angle of the bar (i.e., the m = 2 Fourier component of the 
disk surface density) remains approximately constant upto 
a certain radius and starts varying beyond that. We mea- 
sure the length of the bar as the radius at which the phase 
angle of the bar starts deviating from the contant value. We 
also note that the corotation resonance (R cr = 0.994, 1.12 
and 1.30-Rd at T = 0.56, 1.1 and 2.1 Gyr respectively) of 
the bar is clearly outside the radius confining most of the 
bulge particles. The ratio of R cr /Rbar lies between 1 — 1.4 
at times mentioned in Fig. [U As shown in Fig. [U more and 
more particles in the classical bulge are trapped at the bar's 
ILR as time progesses. During the slow secular evolution, 
we can affirm that the mechanism acting in our system is 
the transfer of angular momentum through resonances. On 
the other hand, during the rapid dynamical evolution, reso- 
nant as well as stochastic orbits played an important role in 
transferring a significant fraction of the net angular momen- 
tum to the classical bulge. This angular momentum transfer 
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Figure 8. Normalized radial profile for the C14 coefficient for the 
classical bulge alone at three different times. 



and the subsequent change of the orbital structure of the 
classical bulge are indeed responsible for the transformation 
of the classical bulge as described below. 



6 EVOLUTION OF THE CLASSICAL BULGE - 
STRUCTURE AND KINEMATICS 

From Sections [2] and [5] we learn that a classical bulge can 
absorb a non-negligible fraction of the total angular momen- 
tum emitted by the bar through resonant interaction. The 
angular momentum gained by the bulge (being a smaller 
mass object than the dark matter halo) has a profound ef- 
fect on its structure, kinematics and dynamics. The result of 
the bar-bulge interaction in our simulation is the transfor- 
mation of an initially non-rotating low mass classical bulge 
into a highly rotating triaxial one. Below we describe, in con- 
siderable detail, various diagnostics which show that this is 
indeed true. 



6.1 Surface brightness 

In the left panels of Fig. [7] we show the surface density 
maps for the classical bulge (viewed edge-on) at four dif- 
ferent epochs during the evolution. The classical bulge is 
shown edge-on (i — 90°) such that the major axis is along 
X-axis and the minor axis along Z-axis. Initially the bulge 
is isotropic and flattened by the strong gravity of the disk 
potential. At later phases of evolution, the inner regions of 
the bulge become rounder and the outer parts become disky. 

In order to understand the structure of the classi- 
cal bulge more quantitatively, we have also performed an 
isophotal analysis using the IRAF ellipse task on a set 
of edge-on images of the bulge including the ones pre- 
sented in Fig. [7] and compute the fourth-order Fourier co- 
sine coefficient ci4 /a normalized to the semi-major axis a 
at which the ellipse was fit. Fig. [S] shows the normalized 
a4 prodiles at three different epochs in the simulation. The 
value s of 04 /a determine the degree of boxiness or diski- 
ness (|Nieto fc Benderll 19891 ). with a^/a < denoting a boxy 
isophote, 04/11 > a disky isophote, and 04 /a ~ means 
elliptical or round isophote. The inner region (X/Rd ^ 0.2) 
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Figure 9. Radial variation of the bulge rotational velocity nor- 
malized to the average velocity dispersion in the central region, 
for five snapshots from T=0 to T=2.1 Gyr. The long tick mark 
on the x-axis denotes the initial value of Ri,, see Table [l] 



of the classical bulge becomes mildly boxy at T = 0.56 Gyr 
and the boxiness increases at 1.1 Gyr as can be seen from 
Fig. and Fig. [7] On the other hand, at T = 2.1 Gyr, 
04 1 a > in the outer parts of the bulge indicating disky 
isophotes. Recall that the disk has a boxy bulge formed as a 
result of the bar buckling instability as shown in Fig. [4] In 
order to compare properties of the classical bulge and boxy 
bulge, we have measured the 04/0 parameter of the boxy 
bulge in a similar fashion as outlined above. It is found that 
at T = 1.1 Gyr, 04/0 of the boxy bulge is negative inside 
X/Rd < 1.0 and its maximum value is about twice that of 
the classical bulge. From the minor axis density profiles cal- 
culated separately for the classical bulge and the boxy bulge 
region, we find that initially the classical bulge extends fur- 
ther above the disk midplane, and its central surface density 
is 2.4 times lower than that of the disk. As the disk goes 
through the buckling instability, the particles settle into the 
3D boxy bulge. We find that at T = 1.1 Gyr the boxy bulge 
is more concentrated toward the disk midplane (z = 0) and 
its midplane surface density is ~ 3 times higher than that 
of the classical bulge. At this time, the density of classical 
bulge above z = 0.177?^ is higher than that of the boxy 
bulge, and the classical bulge extends further. However, at 
T — 2.1 Gyr, the boxy bulge dominates over the classical 
bulge for z ^ 0.39Rd and above this height, their density 
profiles are comparable. 



6.2 Rotational properties 

The influx of angular momentum to the initially non- 
rotating bulge enforces the bulge particles to have a net 
rotational motion. In Fig. [9] we show the radial profiles of 
the rotational velocity normalized to the average velocity 
dispersion in the central region (R ^ R1/2) °f ^ ne classi- 
cal bulge at different epochs during the evolution. The ro- 
tational velocity profiles remain nearly unchanged at later 
stages of evolution when the rate of angular momentum gain 
by the bulge also nearly saturates as can be seen from Fig. [5] 
and Fig. [6] 
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Figure 10. Parallel minor axis velocity profiles of the bulge at 
T=0.56 Gyr (upper panel) and T=2.1 Gyr (lower panel). The 
upper panel shows no clear signature of cylindrical rotation. But 
at later stages of the evolution, cylindrical rotation develops in 
the inner region of the bulge, as indicated by the parallel shapes 
of the velocity profiles in the lower panel. 
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Figure 11. Major axis velocity profiles of the bulge at T=0.56 
Gyr (upper panel) and T=2.1 Gyr (lower panel). The upper panel 
shows no clear signature of cylindrical rotation. But at later stages 
of the evolution, cylindrical rotation develops in the inner region 
of the bulge. The long tick mark on the x-axis denotes the initial 
value of Rf, , see Table [T] 



To illustrate the evolution further, we present four ve- 
locity maps of the classical bulge on the right panels of 
Fig. [7] During the initial phases of secular evolution, the 
angular momentum gained by the bulge particles is pri- 
marily converted into streaming motion and the classical 
bulge starts rotating around the z-axis, with gradients in the 
streaming velocity both along the radial and vertical direc- 
tions. Note that in the barred potential, the classical bulge is 
no longer axisymmetric and its inner regions becomes mod- 
erately boxy (see Sections 16. II and 16. 3|l . This could be a sig- 
nature of a thick bar formed inside the classical bulge (see 
Section [6.3(1 . Indeed as time progresses, mild signatures of 
cylindrical rotation emerge in the inner regions of the bulge 
and gradually become prominent (see Fig. [7} . 

To have a clearer picture of the velocity structure, we 
show parallel minor-axis (Fig. I10|l and major-axis (Fig. 
velocity profiles of the classical bulge at two different epochs 
T — 0.56 (upper panels) and T = 2.1 Gyr (bottom panels). 
The minor-axis velocity profiles are drawn at two different 
radii (X/Rd — 0.31 and X/Rd = 0.71) on either side of 
the bulge center. At T = 0.56 Gyr, the minor axis rotation 
velocity decreases along the vertical direction (^- < 0) in- 
dicating clearly a non-cylindrical rotation throughout the 
bulge (see the upper panel of Fig. I10p . Note, the velocity 
profiles in the outer parts are asymmetric which is probably 
influenced by the on-going buckling instability of the bar. 



The major-axis profiles are taken at four different slits (the 
slit positions are indicated in Fig. Hip parallel to the major 
axis of the classical bulge. The major-axis velocity profiles 
in the upper panel of Fig. 1111 also indicate non-cylindrical 
rotation throughout the classical bulge. 

At later times, the inner regions of the classical bulge 
have developed cylindrical rotation. However, the gradient 
of this cylindrical rotation in the classical bulge is shallower 
than that in the boxy bulge. The minor-axis velocity profiles 
in the bottom panel of Fig. 1101 show clear indication for that 
in the inner regions. The same is evident from the bottom 
panel of Fig. 1111 The major axis velocity profiles at T = 2.1 
Gyr clearly demonstrate that the inner regions (X/Rd < 
0.4) of the classical bulge rotate cylindrically while the outer 
regions beyond about twice the half-mass radius (2 X R\i 2 ~ 
0.45-Rtj) still maintain differential rotation both along the 
radial and vertical directions (-j^- < 0). So in the later stages 
of the secular evolution, the initially non-rotating classical 
bulge has developed a mixed rotational state with the inner 
region rotating cylindrically while the outer region rotates 
differentially in z. 
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6.3 The classical bulge-bar 

In order to achieve a deeper understanding of the complex 
non-linear dynamical interplay of the bar and the bulge, 
we investigate the three dimensional structure of the classi- 
cal bulge using spherical harmonics analysis. In particular, 
we are looking for non-axisymmetric modes in the classical 
bulge which could have been influential for producing some 
of the complex structure and kinematics as discussed in Sec- 
tions l6.ll and l6.2l The outcome of the bar-bulge interaction is 
not only the transfer of angular momentum between the two 
and changes in kinematics thereby, but a structural transfor- 
mation of the classical bulge, a prediction of which is proba- 
bly beyond the scope of the analytic /semi-analytic theories 
|Lvnden-Bell fc Kalnai^ll972l ; iTremaine fc Weinberg||l984r i 
briefly outlined in Section [2] From our analysis, it becomes 
clear that the interaction of a bar and a small classical bulge 
is more vigorous than that between the bar and the massive 
dark halo as discussed in section [5] The primary reason be- 
ing the smaller mass and size of the bulge compared to the 
dark matter halo. 

To analyze the structural components developed in the 
small classical bulge after the evolution, we expand the full 
three dimensional bulge density distribution (p) in terms of 
spherical harmonics: 



1=0 m=-l 



Pi m (r)YT(jS,4>), 



(5) 



where r, 0, <j> are the usual spherical coordinates, and the Y t m 
are the spherical harmonics. p; m denotes the radial density 
function. We bin the bulge particles into spherical shells and 
compute Bi m as function of the bin radius (rt), as follows: 



Bi m (r k ) = Ni m m b Pi m (cos 9j 



(6) 



where Ni m — ((2? + l)/27r) x (l — m)\/(l + m)\, mi, is the mass 
of each bulge particle, P ; m are the Associated Legendre poly- 
nomials. The function Bi m is directly related to the mass of 
each bin and thereby to p; m via the bin radius (r^). Then 
using the above formula (Eq. [6|, we can derive the radial 
variation of the amplitude of a particular I, m mode in the 
bulge as Si m (r k ) = \/5ftB( m 2 + QBi m 2 . The corresponding 
phase angle <j>j can be used to derive the pattern speed of 
the I, m mode. 

In Fig. 1121 we show the time evolution of the amplitude 
of I — 2, m = 2 mode. The classical bulge-bar (hereafter, 
denoted as CIBb) is weaker than the disk bar (see Fig. [3] for 
the bar amplitude and pattern speed) but rotates nearly in 
phase with it. By analyzing the radial variation of the phase 
angles, we conclude that the physical size of the CIBb is 
much smaller compared to the disk bar. Initially, the CIBb 
and disk bar are not in phase, the CIBb seems to be lagging 
behind the disk bar by about 2° — 4° in angle. But soon, 
they start rotating in-phase with each other. After about 1 
Gyr, the pattern speed of the CIBb is also nearly the same as 
that of the disk bar (see Fig. I12|) . A convenient way of view- 
ing the dynamics of the CIBb, is to think of it initially as a 
driven oscillation phenomenon where the disk bar is acting 
as a driver and the bar-like structure in the classical bulge 
is its forced response. Later bulge particles are trapped by 
the 2:1 resonance; i.e., both components populate the or- 
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Figure 12. The strength of the classical bulge-bar (1=2, m=2 
mode) and its pattern speed evolution. 



bits in their jointly rotating pot ential. It has been shown 
in previous simula t ions e . g.. by | Hollev-Bockelmann et al.l 
|2005l ). ICoh'n et all (|200rJ ). lAthanassoulal (|200^ that a bar- 
like structure also forms in the inner regions of the dark 
matter halo as a result of its interaction with the bar in 
the disk. These studies have shown that such a bar in the 
halo is rather weak and nearly corotates with the bar in 
the disk. It turns out that some of the characteristics of 
the CIBb are quite similar to that of the halo-bar. How- 
ever, with the classical bulge being much less massive than 
the halo, the dynamical impact of the bar-like structure is 
much more pronounced in the classical bulge as we have 
already demonstrated above. Beside the transfer of energy 
and angular momentum between the disk-bar and the clas- 
sical bulge, the stars in the classical bulge are also being 
heated during the evolution and hence the inner bulge re- 
gion becomes moderately thicker and rounder (see Fig. [7]). 
We have checked that the slow variation in the ellipticity 
of the classical bulge is consistent with the variation in the 
kinetic energy tensor in accordance with the tensor virial 
theorem. A more detailed picture of the dynamics of the 
bulge hosting a bar and its observational properties will be 
presented in a future paper. 



6.4 Triaxiality and anisotropy 

From the misalignment of the photometric major axis of the 
disk and the bulge and the isophotal twists, it is inferred 
that many o f the bulges in spiral galaxies are indeed tr i axial 
dStarklll977l: ICerhard et a.1 1 Il989l : iBertola et al]ll99ll : lAnnl 
ll995l : lMendez-Abreu et al.ll201fj|) . Here, we show the evolu- 
tion of the triaxiality and velocity anisotropy in the low mass 
classical bulge in our simulation. The global parameter for 
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Figure 13. Time evolution of the shape of the classical bulge. Ini- 
tially the bulge is flattened by the strong disk potential and hence 
oblate. At later phases during the secular evolution it becomes 
triaxial. The solid diagonal line denotes a prolate configuration. 
The red open circles are the measured values of the axes ratios 
at T = 0, 0.25, 0.393, 0.56, 1.1, 1.52, 1.77, 2.1 Gyr. b/a = 0.99 and 
0.877 at T = and 2.1 Gyr respectively. 



I .5 



0.5 



• • . a a 5 . 

. □ . D 



A g 
O 



0.5 



R/R d 



T)=0.5 

-2S 7) £-0.1 



Figure 15. Radial variation of velocity dispersion ratios com- 
puted from bulge particles trapped at ILR (77 = 0.5) and from 
a group of resonant and non-resonant particles in the frequency 
range —2.0 ^5 r\ —0.1 in the classical bulge. Orbits of these 
particles were computed at 1.1 Gyr. 
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Figure 14. Radial variation of the anisotropy parameter rz for 
the classical bulge. Beyond about 0.56 Gyr, the anisotropy pro- 
files remains nearly unchanged. 



the bulg e triaxiality Tj, can be computed using the following 
relation (|Franx et al. I llQQll : IJesseit et al.|[200a ): 



T b = 



1 ~ {b/a? 
1 - {c/af 



(7) 



where a, b, and c are the semi-axes defining the shape of the 
classical bulge (see Section 16.51 for the measurement of the 
axis ratios), a — b > c denotes an oblate configuration i.e. 
Tt — 0, and b = c < a is a prolate figure corresponding to 
Tb = l. a / 6 / c defines a triaxial configuration with peak 
values reaching Tt = 0.5. In Fig. 1 131 we show the evolution of 
the shape of the classical bulge. Initially the bulge is oblate 
(Tb ~ 0, see Fig. I13[) ; thereafter it evolves as a result of the 
angular momentum gain and change in the orbital structure. 
During the period of 0.17 — 0.56 Gyr (roughly the dynam- 
ical phase) a considerable fraction of angular momentum is 
gained at resonances 77 = — 1 corresponding to the Lagrange 
point orbits and r\ — 0.5 (ILR) (see Fig. [BJ. This angular 
momentum gain cause the bulge particles to move outwards 



and thereby producing a disky structure (see Fig. [7J) in the 
outer parts of the bulge. In this period, essentially only b/a 
changes while c/a remains nearly constant. 

Beyond ~ 0.56 Gyr, the CIBb forms in the bulge causing 
substantial changes in the bulge structure. A s mentioned in 
section ROl the CIBb heats |Saha et alj|20ich the bulge stars 
mainly in the central region and makes it thicker. We think 
that this heating due to CIBb is primarily responsible for 
subsequent changes in the c/a. During the period from 1.5 — 
2.1 Gyr, b/a changes more than c/a. At T = 2.1 Gyr, the 
small classical bulge is triaxial with Tt = 0.48. This suggests 
that more generally, fast rotating and triaxial bulges could 
have developed through the interaction of a strong bar and 
a small classical bulge in galaxies with low B/D ratio. A 
more comprehensive analysis focusing on the role of the most 
important parameters such as the bulge mass and size will 
be presented in a future paper. 

As shown in Section |3j the initial velocity distribution 
in the classical bulge in our simulation is isotropic repre- 
sented by a King model. As a result of the angular momen- 
tum influx and the readjustment of the orbits, the veloc- 
ity structure changes during the evolution. We measure the 
deviation from isotropy in the velocity distribution by the 
anisotropy parameters defined as /3 rz = 1 — (a z /a r ) 2 and 
Prtp = 1 — (ov/ov) 2 , where oy, ct^ and a z are the veloc- 
ity dispersions in the radial, azimuthal and vertical direc- 
tion. Then /3 rz > denotes radial anisotropy and p rtp < 
means tangential anisotropy. In Fig. 1141 we show the pro- 
files of the radial anisotropy at four different epochs cal- 
culated from all the particles in the classical bulge. We see 
that the classical bulge already becomes radially anisotropic 
at T = 0.56 Gyr. To understand the source of radial 
anisotropy, we have studied orbits in the classical bulge as 
clarified in Fig. [6] Fig. 1151 shows radial variation of o" 2 /o- 2 
and o^/cr 2 computed from bulge particles that are trapped 
at ILR (77 = 0.5) and from a group of resonant and non- 
resonant particles in the frequency range —2.0 ^ 77 ^ —0.1 
separately. The particles that are at ILR resonance (2:1) fol- 
low xl-type orbits and they produce the radial anisotropy in 
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Figure 16. Vm/cr — e relation for the small classical bulge alone. 
Initially the bulge is non-rotating but then it acquires angular 
momentum emitted by the bar and evolves into a fast rotating 
triaxial bulge. Each open circle represents an epoch in the simu- 
lation beginning at T = (the bottom most point) . Subsequent 
circles are drwan at 0.25, 0.393, 0.56, 1.1, 1.52, 1.77, 2.1 Gyr. The 
pink solid line is the reference isotropic rotator model. 



the inner region of the classical bulge which host the ClBb. 
Whereas in the outer region, definite contributions to the 
radial anisotropy comes from the bulge particles in the fre- 
quency range —2.0 ^ r\ ^ —0.1. The particles that are at 
resonance with the bar e.g., at rj — —2,-1 follow regular 
orbits while the non-resonant particles are in stochastic or- 
bits. 



6.5 Vm/cr — e relation 

To quantify the degree of ordered motion in bulges and el- 
lipticals and illuminate the difference between the two types 
of stellar systems, the Vm/cr— e diagram relating the ratio of 
rotational to ra ndom motions and the observed ellipticity (e) 
was introduced |lllingworthll 19771 ) . It was shown that bulges 
are, in g eneral, fast rotators compared to bright elliptical 
galaxies (iKormendv fc Ilhngwortl]|l982l; iDavies i et al.l l l9S3l ; 



ICappellari et al.ll2007l ; iMorelli et al.ll2008l ). Here, we focus 
on the relation between the shape and the kinematics of the 
simulated low mass classical bulge that has been subject to 
the secular evolution driven by a strong bar. We show explic- 
itly the evolutionary track of this particular classical bulge 
in the V m /a — e diagram below (Fig. 116}) . 

In observations, it is rather difficult to have an accurate 
measurement of the bulge rotation velocity due to possible 
disk contamination. On the other hand, in simulations, it 
is rather straightforward to compute the velocity profile for 
the classical bulge alone because it is possible to filter out 
the disk and halo components of the model galaxy. We de- 
termine, Vm. as the maximum of the azimuthally averaged 
rotational velocity of the bulge particles measured in the 
equatorial plane of the bulge; a is the mean velocity disper- 
sion in the central region (calculated at r ~ 0.5 x R1/2) of 
the bulge. 

For the bulge ellipticity, we use c/a for edge-on view. 
Measuring the ellipticity for a bulge is a bit tricky because 
there can be strong radial variation in the ellipticity pro- 



file e(r). Triaxiality could add another degree of complexity 
to such a measurement. Below we describe how the bulge 
ellipticities are measured. 

In order to determine the intrinsic ellipticity, we first 
compute the moment-of-inertia tensor of the three dimen- 
sional mass distribution of the classical bulge and diagonal- 
ize it to obtain the principal moments and three orthogonal 
eigenvectors. The principal moments of inertia determine the 
intrinsic axis ratios of the inertia ellipsoid and the eigenvec- 
tors determine the orientation of the ellipsoid with respect 
to the co-ordinate space. Using the eigenvalues and eigen- 
vectors, we determine the two axis ratios namely b/a and 
c/a, where a > b > c are the three semi-axes of the inertia 
ellipsoid. We have done this at the bulge half-mass radius 
containing 50% of the total bulge particles and at the radius 
containing about 90% of the total bulge particles. The el- 
lipticity measurements for the classical bulge are nearly the 
same in both cases. In the following, we use the ellipticity 
at radii enclosing ~ 90% of the total bulge particles and use 
the definition of the ellipticity of the bulge as e = 1 — c/a 
when viewed edge-on. We have also performed isophote anal- 
ysis using the ellipse-fitting routine from IRAF on a set of 
suitably rotated and inclined edge-on images (Fig. [7} of the 
classical bulge at different epochs during the evolution. We 
find a good agreement between the two different types of 
measurements of the bulge ellipticity. 

In Fig. 1161 we show the Vm/a and e values for the small 
classical bulge during the secular evolution. Each point in 
this diagram corresponds to a particular epoch during its 
evolution and when connected together they form its evo- 
lutionary track. This shows that the small classical bulge 
rotates significantly faster in the latter stages of evolution 
compared to the obl ate isotropic rotator model, wh ich can 
be approximated by (|Binnevlll978l ; lKormendvlll982r i: 



V m /o £ 



(8) 



The interpretation of Fig.[l6]is complicated by the fact that 
the classical bulge is not an isolated stellar system, but inter- 
acts dynamically with the bar within a disk galaxy. The en- 
tire period of evolution of the classical bulge can be broadly 
divided into two parts: one before the formation of the ClBb 
(~ 0.56 Gyr) and the second after its formation. Before the 
formation of the ClBb, the bulge stars in the outer region 
gain a significant fraction of angular momentum emitted by 
the bar and move mainly outwards in radius. As a result 
of this, the values of Vm/cr increase till ~ 0.56 Gyr while 
the axis ratio c/a remains unchanged. In the second half of 
the evolution, the bulge stars are heated due to the ClBb 
by a factor of ~ 1.5 in velocity dispersion and this makes 
the inner region rounder. So c/a ratio increases, making the 
ellipticity decrease considerably. The near saturation in the 
Vm/cr towards the end of the simulation is connected with 
the fact that the rate of angular momentum gain by the 
bulge nearly saturates at these epochs (see Fig. [5] Fig. [5} . A 
more detailed analysis on how the spinning up of the clas- 
sical bulge depends on the various parameters of the bulge 
(bulgc-to-disk mass ratio, size of the bulge, its central veloc- 
ity dispersion) and disk (Toomre Q, bar strength, bar size) 
will be presented in a future paper. 
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7 DISCUSSION 

The primary goal of this paper has been to describe a generic 
mechanism, the transfer of angular momentum from a bar 
to an embedded classical bulge. We have shown that this 
mechanism is important for understanding the rotational 
motion of low mass classical bulges in the central regions 
of barred galaxies. The growth rate and the strength of the 
bar are important factors for the mechanism to work effi- 
ciently. One very interesting outcome of this process is the 
cylindrical rotation in the small classical bulge in the model 
studied here. Some possible observational implications and 
other issues are addressed below. 



classical bulge. The upper limit on the total mass of the 
bulge (boxy bulge -fclassical bulge) in our model, including 
the remaining disk component in the boxy bulge region, is 
~ 1.46 x 1O 1O M . Of this, 0.29 x 1O 1O M is in the classical 
bulge, and ~ 1.17 x 1O 1O M0 in the boxy bulge and the cen- 
tral disk. Since the classical bulge extends further above the 
galactic plane than the boxy bulge, the metallicity composi- 
tion of the composite system would change with height. We 
plan to investigate this further and use our model to look for 
signatures of the classical bulge from the metallicity distri- 
bution in order to compare with observations of the Milky 
Way. 



7.1 Structural impact on the bulge 

During the secular evolution, a flat bar transforms into a 
boxy bulge. While the bar grows, buckles and evolves, a 
fraction of the angular momentum emitted by the bar is 
absorbed mainly in the outer parts of the embedded classi- 
cal bulge in the galaxy. As a result, streaming motions are 
induced in the classical bulge (Fig. [9}, and the orbital struc- 
ture changes, causing velocity anisotropics. After about half 
a Gyr, the CIBb forms in the inner regions of the bulge. The 
CIBb introduces a pattern rotation in the classical bulge 
which transforms into a rotating triaxial object. A compar- 
ison of Fig.[3]and Fig. H2l shows that beyond ~ 0.8 Gyr the 
CIBb essentially corotates with the boxy bulge formed out 
of the disk bar. Towards the end of the evolution, the model 
galaxy has a composite bulge: a superposition of the boxy 
bulge and a rotating triaxial classical bulge. At late times, 
the inner regions of the classical bulge become rounder as 
it is evident from the surface density maps (Fig. [7]) for the 
classical bulge particles alone. We note that the stars in the 
classical bulge are heated by a factor of ~ 1.5 in velocity 
dispersion within 2 Gyrs. The CIBb maybe responsible for 
heating the bu lge stars in a sim ilar fashion as the bar heats 
the disk stars (|Saha et alj|2010h . 



7.2 Boxy bulge and net cylindrical rotation 

It is widely accepted that cylindrical rotation is a character- 
istic feature of the kinematics of a boxy bulge formed out 
of disk material, via the vertical buckling instability of the 
bar. Thus the presence of cylindrical rotation in the central 
regions of a galaxy may lead one to infer the presence of a 
boxy bulge that originated from the disk, without any need 
for a classical bulge in this galaxy. 

The work presented in this paper adds a new aspect to 
this simple picture. The cylindrical rotation could also in- 
clude the stars of the classical bulge whose rotational prop- 
erties have been modified by the interaction with the bar. 
One would measure the net cylindrical rotation of the stars 
in the combined boxy bulge and classical bulge. In absence 
of strong photometric evidence, other information such as 
from stellar populations and metallicity gradients would be 
needed to determine the presence of a small classical bulge. 

A lthough the pure kinematic modelling of the BRAVA 
data |Shen et al.l l2010h suggests only an upper limit on 
the mass of a classical bulge in the Milky Way, the mea- 
surem ents of the metalli city gradient above the galactic 
plane (|Zoccali et al.l [20081 ) may indicate the presence of a 



7.3 Observing secular evolution through a 
classical bulge 

Previous studies mainly based on iV-bod y simulations 
have focused on the b a r-halo i nteraction (jAthanassoulal 
20021; IWeinberg fc Katzl l2007al ; IVilla- Vargas et al.l |2009| ; 



Dubinski et al.N2009l 'l and shown that a significant amount 



of angular momentum emitted by the bar is absorbed 
in the dark matter halo. The angular momentum gained 
by the halo changes the internal structure of the dark 
matter halo. Some authors have utilized this mecha- 
nism of ang ular momentum exchange to resolve the cusp- 
core issue llEl-Zant et al.l 1200 ll ; IWeinberg fc Katzl [2002; 
ISellwood| [2008) in galaxies while others have focused on the 
halo-bar (iHollev-Bockelmann et aI1l2005l ; ICoh'n et al.ll2006l ; 
lAthanassoulal 20071 ). However, direct observational evidence 
for the halo-bar, and hence direct observational verification 
of the ongoing secular evolution and angular momentum 
transfer can not be obtained unless dark matter is detected. 
Unlike the dark matter, it is possible to observe galactic 
bulge stars in detail, both the kinematics and stellar popu- 
lation parameters. It is thus possible to verify observation- 
ally the bar-bulge interaction and the resulting dynamical 
properties of the CIBb. Observational evidence for the CIBb 
could be a direct confirmation of the angular momentum 
transfer and secular evolution in the galaxy. 



8 CONCLUSIONS 

The secular processes driven by the bar not only restruc- 
ture the disk, but also the other components in the galaxy. 
Since the classical bulge is less massive here compared with 
the surrounding dark matter halo, the angular momentum 
gained by the classical bulge has a more significant effect 
on its evolution. The present work has shown in consider- 
able detail that the unavoidable gravitational interaction be- 
tween these two components can have profound implications 
for the structure of a low mass classical bulge, as highlighted 
below. 

1. We have established that the main mechanism of an- 
gular momentum transport operating between the bar and 
the classical bulge is through resonances. The bulge parti- 
cles gain angular momentum emitted by the bar through the 
bar's ILR (?y = 0.5), and other resonances with r\ = —1, —2 
and also through non-resonant orbits with r\ G (—0.4,0.) 
during the dynamical phase when bar growth is rapid. Ap- 
proximately 3/4 of the net angular momentum is gained 
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by the classical bulge during the dynamical phase where 
stochastic orbits contribute ~ 30%. 

2. The angular momentum gained by the initially non- 
rotating classical bulge sets the bulge particles in rotational 
motion. The radial gradient in the rotational motion in the 
classical bulge is lower than in the boxy bulge. As time pro- 
gresses the rotational velocity increases and nearly saturates 
at about 1 Gyr. At around this time, the inner regions of the 
classical bulge develop cylindrical rotation while the outer 
parts are still in differential rotation. 

3. As a result of the angular momentum transfer, 
some of the bulge orbits are trapped by the rotating bar 
potential, a CIBb forms which essentially cororates with 
the bar, and the classical bulge transforms into a triaxial, 
anisotropic object where trapped orbits contributes to the 
radial anisotropy. 

4. Towards the end of the secular evolution, the model 
galaxy has a composite bulge which is a superposition of the 
boxy bulge formed out of the disk material and the rotating 
triaxial, low mass classical bulge. The stars in the compos- 
ite bulge rotates cylindrically. From an observational per- 
spective, one would need other tracers such as metallicity 
gradient, stellar population parameters along with the kine- 
matics to reliably determine the presence of such a low mass 
classical bulge embedded in the boxy bulge. 
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